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ABSTRACT 

The objective of this project was to develop an interac- 
tive computer program providing flexibility in the design 
and analysis of marine propulsion gears. The program/ 
Reduction Gear Analysis and Design (R2GAD) , will handle con- 
vential parallel axis and simple epicyclic reduction gears. 
It is capable of generating preliminary designs of new gear 
sets or providing analyses of existing or proposed gear 
sets. Program development, organization, and operation are 
discussed. 
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I 



INTRODUCTION 



In the conceptual stage of ship design, many parameters 
and options are considered. This is especially true with 
respect to the propulsion plant. Changes in hull design, 
displacement, and numbers of propellers all affect changes 
in the requirements for the propulsion plant. There are 
also various options under consideration in the propulsion 
plant: turbines or internal combustion engines, the number 
of enqines, the auxiliaries required for support, etc. All 
of these will affect the initial design of the reduction 
gears. It is, therefore, important to be able to produce 
preliminary designs of reduction gears for the options under 
considerat ion . 

Preliminary designs provide useful information on feasi- 
ble arrangements and size without going into the specific 
design details dependent upon manufacturing. Since any num- 
ber of preliminary designs may be required due to perturba- 
tions discussed above, it is important to be able to 
automate the design process. An interactive computer pro- 
gram providing various options would free the engineer from 
tedious, time consuming, and often error prone number manip- 
ulation and allow him to produce multiple designs for con- 
sideration. If would also provide a quick means of checking 
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the affects of various parameters in addition to the ability 
to analyze proposed designs or configurations. 

Reduction Gear Analysis and Design (R2GAD) was developed 
to fill this need. It is an interactive computer program 
offering close user control through numerous options. Being 
interactive, it provides a rapid means of designing or ana- 
lyzing a gear set, thereby reducing the turn-around time 
inherent in the use of batch systems. The program was kept 
modularized and well documented for ease of maintenance and 
modification. The modularized construction also provides an 
additional benefit of being able to use this program on 
smaller computers by using an overlay scheme. 



II. PROGRAM CAPABILITIES 



A. SCOPE 

REGAD was written to provide preliminary designs or 
analyses of marine propulsion reduction gears. It is inca- 
pable of providing detailed designs or performing detailed 
analyses since specifics of manufacturing are not reguired 
for input. The program does not consider shafting, bear- 
ings, lubrication, couplings, casings, or other auxiliaries. 
It will provide sizing information in the form of pitch 
diameters, effective facevidths, gear ratios, and numbers of 
teeth per gear. In addition, the program will provide esti- 
mates of loadings and stress levels. Estimated weight and 
dimensions of the gear set are also provided. 

All computations are based primarily on the American 
Gear Manufacturers Association's standards [Ref. 1, 2, 3] 
using appropriate constants for marine propulsion gears 
[Ref. 4, 5]. As an option in the program, these constants 
can be replaced by the user to enable him to investigate 
other applications such as reduction gears for ships service 
or emergency generators. 
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3. LIMITATIONS AND OPTIONS 



Program application is limited to marine reduction gears 
with a maximum of three reduction stages. Conventional par- 
allel axis and simple epicyclic arrangements with helical 
gears are possible. When dealing with epicyclics, it is 
assumed that load sharing of the planets is achieved and 
that the ring gear is suitably flexible. Efficiencies of 
the gear sets are not provided since power losses are not 
computed. While estimates of bending and contact stresses 
are provided, scoring can not be estimated since lubrication 
is not considered. REGAD does not require the K-factor as 
input as in previous programs since hardness ranges for the 
pinions and gears are required. However, the K-factors are 
computed and displayed for reference purposes. The weight 
estimates are based on actual designs and do not include 
turning gears, attached lubrication oil pumps, or other 
auxiliaries. 

The following is a list of major options provided by the 
program : 

(1) brief, on-line program description 

(2) choice of design or analysis 

(3) listing of preprogrammed constants and an ability to 
change selected constants of the user's choice 

(4) choice of single, double, or triple reductions 

(5) choice of single or double helical gears 

(6) choice of six hardness ranges for gears and pinions 
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(7) conventional parallel axis arrangements (see Figures 
1 , 2 , and 3) 

(a) one or two power inputs 

(b) single power path (articulated) or dual power 
paths (locked train) 

(8) simple epicyclic arrangements (see Figure 4) 

(a) choice of planetary or star arrangements 

(b) single power input 

(c) choice of three, four, or five planet/star gears. 
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III. PROGRAM ORGANIZATION AND OPERATION 



A. REGAD FLOWPATHS 

As stated previously, the program was designed in modu- 
lar form with each module consisting of a number of subpro- 
grams. These modules are just conceptual groupings of 
associated subprograms, and are not related to actual pro- 
gram implementation on any specific computer. Figure 5 
shows the basic flow paths of the program. Module One is 
for program initialization and problem set up. Module Two 
performs calculations for conventional parallel axis gear 
sets, while Module Three handles epicyclic gear sets. Mod- 
ule Four is a grouping of all the computaxional subprograms 
required by the other modules. 

B. MODULE IDENTIFICATION AND DESCRIPTION 

This section provides a brief description of each sub- 
program in each module. 

1 • Module One : Initialization and Set-up 

Module One contains the subprograms necessary for 
initialization, execution, and initial data entry. It is, 
basically, the control module for the program. The follow- 
ing is a grouping of the subprograms in Module One. 
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a. P EG AD 

HEGAD is the main program. It provides the 
options for either design or analysis and either parallel 
axis or epicyclic arrangements and controls the flow to the 
proper module. It then calls the required subprograms for 
execution of Modules One, and Two or Three. 

b. BLOCK DATA 

The BLOCK DATA subprogram initializes variables 
in each of the common blocks. 

c. SUBROUTINE DSCRPT 

This subroutine is called by REGAD after an 
affirmative response to a user option to provide a brief 
description cf the REGAD package. It contains an option to 
stop the program if only a program description is desired. 

d. SUBROUTINE INPUT 

All options and initial design parameters are 
entered via this subroutine which is called by REGAD. 

e. SUBROUTINE AGMA 

The constants for marine propulsion gears 
required by various AGMA formulations are initialized in the 
BLOCK DATA subprogram, and can be listed as an option in 
REGAD. REGAD calls this subroutine after an affirmative 
user response to display the preprogrammed values. This 
subroutine then allows the user to selectively change any 
constant desired. 
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2 . 



Module Two : Parallel Axis 



This module contains all the major subprograms 
called by REGAD to provide an initial design or to perform 
an analysis of convential parallel axis reduction gears. 

The following is a grouping of the subprograms in Module 
Two . 

a. SUBROUTINE PRLDES 

This subroutine will produce a design of a par- 
allel axis gear set. All pinion and gear diameters, effec- 
tive facewidths, and gear ratios are computed using a basic 
random search optimization technique to find a feasible 
design by attempting to minimize a function of gear pair 
volume. It should be noted that, while attempting to mini- 
mize gear volume, the design is not necessarily optimized 
for minimum weight. The optimization technique is used here 
only to produce a feasible design in terms of dimension and 
power constraints by minimizing a function of gear pair vol- 
ume. To produce a truly optimized design for minimum 
weight, a full optimization must include many more design 
variables such as helix and pressure angles, pitches, and 
hardnesses in addition to the dimensions. Additional con- 
straints such as stress and unit load levels would need to 
be incorporated. All of this would require a more sophisti- 
cated and efficient optimization technique than is used 
here. 
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b. SUBROUTINE PRLANL 



To analyze a proposed or existing design, REGAD 
will call this subroutine. It will request, as user-sup- 
plied input, the basic information calculated in PRLDES, 
i.e., pitch diameters and effective facewidths. Using this 
information, PRLANL will compute other parameters such as 
gear ratios, power and speed splits, and numbers of teeth 
per gear. 

C. SUBROUTINE PRLRES 

Immediately following a call to PRLDES or 
PRLANL, REGAD will call PRLRES to compute all remaining 
information such as expected loadings and stress levels. 

The user should be aware that the stress levels are computed 
according to AG KA formulations [Ref. 2 and 3] and take into 
account load distribution and overloads. This will produce 
levels that may seem high but are actually closer to actual 
levels to be expected in service, 
d. SUBROUTINE PRLSIZ 

REGAD calls this subroutine after PRLRES to com- 
pute estimates of gear set weight and gearbox dimensions. 
These estimates are determined by empirical relationships 
obtained from a rather limited data base of actual designs. 
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e. SUBROUTINE PRLOUT 

This is the last subroutine called by REGAD in 
the parallel axis path. It provides a detailed output of 
the results obtained from the design or analysis including 
design parameters entered by the user, the dimensions of 
each component, expected loadings and stress levels, and 
configuration information. 

3. Module Three : Epicvclic 

Module Three contains all the major subprograms 
called by REGAD no design or analyze simple epicyclic reduc- 
tion gears. Subroutines EPCDES, EPCANL, EPCRES, EPCSIZ, and 
2PC0UT are all analogous to those in Module Two. They per- 
form the same functions, but for simple epicyclic gears. 
Therefore, individual descriptions will not be repeated 
here . 

4. Module Four : Computational Subprogram Library 

This module is an organizational grouping of all rhe 
subprograms called by those in Modules One, Two, and Three, 
a. Subroutine Subprograms 

The following are the subroutines used: 

(1) GFI - subroutine to compute the AGMA dur- 
ability geometry factor, I 

(2) GFJ - subroutine to compute the AGMA 
strength geometry factor, J. 
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b. Real Function Subprograms 

The following are the real function subprograms 

used: 

(1) ARCCOS - computes the arc cosine of two 
arguments 

(2) ARCSIN - computes the arc sine of two argu- 
ments 

(3) AGMAE1 - uses LaGrangian interpolation of 
Table E- 1 [Ref. 1] to compute the constants 
required for the stress concentration factor 
formu lat ion 

(4) CKDATA - called by SUBROUTINE AGMA to allow 
the user to change the preprogrammed con- 
stants 

(5) POWERS - computes allowable service power 
based on AGMA strength rating [Ref. 3] 

(6) POWERH - computes allowable service power 
based on AGMA durability rating [Ref. 2] 

(7) RTFNDR - a modified version of FUNCTION 
ZEROIN [Ref. 6] used to find a zero of a 
function in a specified interval 

(3) FALFA - the function required by SUBROUTINE 
GFJ and the zero of which is computed in 
FUNCTION RTFNDR 

(9) SHRLD - computes the load sharing ratio, m N 
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(10) THICK - computes tootn thickness at any 

diameter given a known thickness at a dif- 
ferent diameter. 

C. DATA TRANSFER 

All data transfer between subprograms in Modules One, 
Two, and Three is via combinations of seven common blocks. 
Data transfer to and from subprograms in Module Four is via 
argument lists and common blocks as required. The following 
is a list of the common blocks used: 

(1) /AGMA3/ : constants for AGHA strength formulations 

(2) /AGMAH/ : constants for AGMA durability formulations 

(3) /DESDAT/ : design parameters and options 

(4) /DESPRL/ : parallel axis design information 

(5) /RES ?RL/ : parallel axis computational results 

(6) /DESEPC/ : epicyclic design information 

(7) /RESEPC/ : epicyclic computational results. 

The variables in each common block along with their defini- 
tions can be found in Appendix 3. 

D. PROGRAM OPERATION 

REGAD is an interactive program designed to allow the 
user to solve his problem at a terminal. Being interactive, 
the program has many options that control program execution, 
in addition to requests for data necessary for the execution 
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of the program. Each request for information will contain 
the necessary guidelines needed by the the user to respond. 
This may take the form of a mini-table containing informa- 
tion on each option choice, the range of values when a spe- 
cific quantity is requested, or units, where applicable, of 
the requested data. 

All option parameters are integer values and should not 
be entered with a decimal. Option codes entered by the user 
are checked for validity to ensure they fall within the 
allowed range. If two options are offered, enter a 1 or a 
2. Any value entered less than one will automatically 
default to one, and any value greater than two will automat- 
ically default to two. In those cases where there are more 
than two options, the response is checked to see if it falls 
within the allowed range. If it does not, a message alerts 
the user to this fact and allows him to re-enter the correct 
code. Some questions require affimative or negative respon- 
ses. To reply, use a Y for yes or an N for no. Use of 
other values may give undesirable results. 

Every attempt has been made to anticipate possible error 
conditions. If one of these is encountered, a message is 
generated to inform the user. If the error encountered is a 
terminal error, the message will also indicate that the pro- 
gram run was aborted under program control. 
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A detailed development of this package is provided in 
Appendix A where specifics can be found. Appendix 3 pro- 
vides a cross-reference of the variables used in Appendix A 
with those used in the program. It also contains detailed 
information on the common blocks. Sample runs of the pro- 
gram can be found in Appendix C, and a complete listing of 
the program is in Appendix D. 
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IV. CONCLUSIONS AND RECOMMENDATIONS 



Computer aided design (CAD) is an important and useful 
tool for engineers. As computer technology continues to 
expand, CAD will become increasingly available for the prac- 
ticing engineer, allowing him to use his initiative in 
design instead of being a slave to the numbers involved. 
REGAD is such a tool for use in the preliminary design of 
marine reduction gears during the conceptual stages of pro- 
pulsion plant design. 

REGAD could become even more useful if additional 
options are provided. A module to perform sensitivity anal- 
yses of a given design would greatly enhance the use of this 
program. This option would allow the user to start with any 
design and vary a selected variable over a specified range 
to determine its impact on the design. It could also be 
used to "fine tune" a design by modifying selected parame- 
ters to produce the results desired without having to rerun 
the program for each modification. Graphics would add 
another dimension by providing graphical displays of the 
gear arrangements and of certain data such as the results of 
a sensitivity analysis. A module to handle various compos- 
ite designs of parallel axis and epicyclic gears would be an 
important addition. Also, it is recommended that a larger 
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data bass be collected to provide more accurate empirical 
constants for the weight and gearbox size estimates. 
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APPENDIX A 



PROGRAM DEVELOPMENT 



With the exception of several general conversion rela- 
tionships, all computations are accomplished in Modules Two 
and Three with calls to subprograms in Module Four. The 
analytical relationships used in the program will be exam- 
ined, however, most of the relations used can be easily 
found in the literature and in various texts, so background 
developments will not be given. 

I. GENERAL RELATIONSHIPS 

The following relationships are used in Module One and 
in various other subprograms. The transverse diametral 
pitch of any gear is the ratio of its number of teeth to its 
pitch diameter; 

N 

P d = (D 

d 

The normal and transverse diametral pitches are related by; 

P d = PndCOS Ip (2) 

and the pressure angles by; 
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tan 0„ = tan 0, cos ip 



( 3 ) 



Axial pitch is defined as; 



n 

?x = 

Pnd Sin 



TT 

P d tan ip 



( 4 ) 



II . CONVENTIONAL PARALLEL AXIS FORMULATIONS 
Subroutines PRLANL and PRLDES each provide the pitch 
diameters of the pinions and gears, the effective face- 
widths, the stage reduction ratios, the numbers of teeth per 
gear, speed and power splits, and the geometry factors to 
subroutines PRLRES and PRLSIZ to compute all further infor- 
mation. The speed splits are the actual speeds of the indi- 
vidual gears and a power split is the actual power 
transferred by a gear. The strength and durability geometry 
factors are computed in separate subroutines in Module Four 
and will be discussed later. 

A. COMMON RELATIONSHIPS 

Power splits are determined from the configuration. For 
a single power path configuration, the power is transferred 
equally from the pinion to the gear, where in a dual power 
path configuration, the pinion transfers one half its power 
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to each of two gears. These splits are computed exactly 
since losses are neglected. 

Speed splits and stage reduction ratios are based on; 

D n G 

m G = = (5) 

d n p 

Numbers of teeth on each gear are computed from the 
equation below and are rounded to the nearest integer. 

N = d x ? d (6) 

B. DETERMINATION OF DIAMETERS AND FACSWIDTHS 

All diameters and facewidths are entered by the user in 
subroutine PRLANL. Stage gear ratios, power and speed 
splirs and numbers of teeth per gear are computed as dis- 
cussed in the previous section. 

In subroutine PRLDES , the diameters, facewidths, and 
stage gear ratios are determined by using a basic local ran- 
dom search optimization technique to produce a feasible 
design. This algorithm requires an initial design to start. 

The initial design is based on Dudley's [Ref. 7] formu- 
lation for preliminary estimates of gear size; 

31500 Pwr (m G + 1 ) 3 

C 2 F = (7) 

K n p m G 
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c 



d 



(iEg + 1 ) 



(8) 



2 

By substituting equation 8 into equation 7, a formula for 
estimating pinion diameter is obtained; 

126000 Pwr (m G + 1) 

d 3 (9) 

n P K (F/d) m G 

where F/d = 1.0 for single helical gears and F/d = 2.25 for 
double helical gears. The term K is the K-factor which is 
an indication of durability. An expression for estimating K 
is provided by Thoma [Ref. 4]; 

s ac x 10-A* /2.80 \ 

— — ) x (tt) (10> 

where the constants used are the AGMA durability constants. 
The K-factor in equation 10 is for the second reduction. 

For the first reduction, multiply K from equation 10 by 
1.20. The initial estimates for the stage gear ratios are; 
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(1) single reduction m G = M 0 



(3) triple reduction 



(2) double reduction 



+ 3 dual power path 
n>0 2 =Vir o - 1 single power path 
m Gi =VF 0 /m G2 
n> G =VH; 

® G3 ='^To + 3 

VF 0 . 



The initial facewidths used are: 

(1) single helical gears F = d 

(2) double helical gears F = 2.25 d. 

With this initial design as a starting point for the 
random search algorithm, successive designs are determined 
by randomly adding small amounts of between +1.0 and -1.0 to 
the diameters, facewidths, and stage gear ratios. These 
small amounts are scaled to take into account the difference 
in range of values for each variable. This process will 
attempt to find a feasible design in which all specified 
constraints are satisfied. If the initial design violates 
one or more constraints, the design that violates them the 
least in succeeding iterations will be kept until a design 
satisfying all constraints is found. Once a feasible design 
is found, an attempt to improve this design is made by try- 
ing to reduce the size of the gears by minimizing a function 
of gear pair volume; 



Volume = £Lc 2 F = £ £ [ i (m G +1) zd2F] 



( 11 ) 
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The interior summation is over the number of reduction 
stages, and the exterior summation is over the number of 
power inputs. The constraints imposed which determine the 
limits on each of the designs are: 

(1) actual transmitted power is less than or equal to tne 
allowable service power in accordance wirh references 
2 and 3 

(2) maximum gear diameter of 200 inches due to manufactur 
ing limitations 

(3) minimum facewidth greater than four axial pitches to 
ensure proper helical action 

(4) maximum facewidth less than the pinion pitch diameter 
for a single helical gear or 2.25 times the pinion 
pitch diameter for a double helical gear 

(5) pinions and gears in succeeding reduction stages are 
to be larger than those in the previous stage due to 
the greater amounts of torque carried 

(6) in dual power path arrangements, the gear ratio for 
each reduction stage is greater than the preceeding 
stage due to the torque carried. 

The design obtained can than be adjusted by the user as 
desired by changing parameters with the analysis option. 
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III. 



5PICYCLIC FORMULATIONS 



As in Module Two, the pitch diameters, effective face- 
widths, stage reduction ratios, numbers of teeth, speed and 
power splits, and the geometry factors are all entered or 
computed in the analysis or design subroutines (EPCAUL or 
EPCDES) for use in the final computations subroutines 
(EPCRES and EPCSIZ) . Here, the speed splits are the rota- 
tional speeds of the sun and planet gears and of either the 
ring gear or the carrier, depending on the configuration. 
Planetary arrangements have fixed ring gears while star 
arrangements have fixed carriers. Also, the direction of 
rotation must be considered. Star arrangements reverse the 
direction of rotation of the input and the planetary 
arrangements will maintain direction of rotation. Assuming 
equal load sharing of the planets and neglecting losses, 
power splits are straightforward. The input and output pow- 
ers are equal while each planet carries an equal share of 
the total power. Load sharing is an important consideration 
in the design of epicyclic gears, and must be assured in 
marine reduction gears due to the high power levels experi- 
enced. Equal load sharing of the planets can be reasonably 
achieved in several different ways. One method requires the 
sun gear to float, supported only by the planet gears, with 
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a relatively flexible ring gear to allow for inaccuracies in 
the teeth. There are also mechanical devices available to 
assist in achieving an equal division of the load. Experi- 
ence has shown, for marine applications, that three to five 
planets with stage ratios in the range of two to eight work 
best . 



A. COMMON RELATIONSHIPS 

Unlike conventional parallel axis arrangements, there 
are specific numerical rules governing the proper assembly 
and operation of an epicyclic gear set. These involve the 
selection of the numbers of teeth and planets along with 
computing the various speed ratios. Mesh frequencies are 
also configuration dependent as seen in a following section. 

There are basically three relationships that must be 
satisfied to ensure proper assembly and operation. The 
first is a relationship defining the speed ratio of the epi- 
cyciic stage since it is not merely the ratio of numbers of 
teeth or diameters as in a conventional parallel axis gear 
set. The second relationship requires the ring gear diame- 
ter to be equal to the sum of the sun gear diameter and 
twice a planet gear's diameter. This ensures the planets' 
ability to fit between the sun and ring gear. For the final 
relationship, it can be shown geometrically that the sum of 
the numbers of teeth on the sun gear and ring gear must be 
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an integral multiple of the number of planets in the gear 
set to ensure proper alignment and meshing of all teeth. It 
should be noted that these relationships are based on 
equally spaced planets around the sun gear. The above rela- 
tionships are conveniently expressed in terms of numbers of 
teeth on each gear as seen in references 7 and 8. The speed 
ratio for a planetary arrangement is; 

n o 

m G = = ♦ 1 (12) 

Hi N s 

and for a star arrangement; 



-n 0 -N„ 

m G = = (13) 

n j Ng 

where the negative sign indicares the star arrangement's 
reversal of rotational direction of the input. The rota- 
tional speed of the planet gears is required for the design 
of their bearings and can be determined by; 

Nr 

n PLN = n 0 (14) 

^PLN 

where n 0 in each equation above is the speed of the carrier 
for a planetary arrangement or is the ring gear's speed for 
a star arrangement. The assembly and meshing relations in 
terms of tooth numbers are; 
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(15) 



Nr = N s + 2 N pln 



and; 



N R + N s = k NP (16) 

where k is an integer and NP is the number of planets. 

B. DETERMINATION OF DIAMETERS AND FACEW IDTHS 

The random search technique discussed for conventional 
parallel axis gears is used to provide the epicyclic diame- 
ters, facewidths, and stage gear ratios. Equation 9 is used 
to provide an initial estimate of sun gear diameters where m G 
is replaced by the ratio of the planet's pitch diameter to 
the sun's pitch diameter. This value is usually in the 
range of 1.5 to 3; therefore, a random number in this range 
is used to start the problem. Once the sun gear diameter is 
estimated, the other diameters can be found using the rela- 
tionships in equations 12 to 16. The initial estimates for 
the stage gear ratio are the roots of rhe overall ratio cor- 
responding to the number of reduction stages. For example, 
m Gi and m G2 for a double reduction gear set would be the 
square root of the overall ratio. Initial facewidths are 
chosen as before. The initial estimates of the diameters, 
facewidths, and the gear ratios provide a starting point for 
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the random search algorithm discussed previously. Again, 
the method will attempt to improve feasible designs by mini- 
mizing a function of gear volume; 

Volume = (NP- d^ N + d| + djj ) • F (17) 

where the summation is over the number of reduction stages. 
The constraints imposed are similar to those for the paral- 
lel axis gears; 

(1) actual transmitted power is less than or equal to the 
allowable service power in accordance with references 
2 and 3 

(2) maximum ring gear diameter of 150 inches due to manu- 
facturing limitations 

(3) minimum facewidth greater than four axial pitches to 
ensure proper helical action 

(4) maximum facewidth less than the sun's pitch diameter 
for a single helical gear or 2.25 times the sun's 
pitch diameter for a double helical gear 

(5) planet gears are to be larger than sun gears due to 
the greater amounts of torque carried 

(6) stage gear ratios are to be between 2 and 8 for each 
reduction stage. 

As before, once a design is obtained, the user can utilize 
the analysis option to obtain the desired results. 
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IV. 



COMPUTATIONAL 3 ESDLTS AND DESIGN INFORMATION 



Once the geometry is determined in the analysis or 
design subroutines, the computational results subroutine 
(PRLRES or EPCRSS) and the size estimates subroutine (PRLSIZ 
or EPCSIZ) are called to provide design information concern- 
ing tooth loads, stresses, and other configuration, geom- 
etry, and size information. This section describes the 
formulations used. 

The facewidth to diameter ratio is computed using the 
effective facewidth and the pitch diameter of the pinion for 
parallel axis gears or the sun gear for epicyclics. Center 
distance is taken as the average of the pinion and gear 
pitch diameters. A center distance is computed for epicycl- 
ics by finding the average of the sun and a planet gears' 
pitch diameters. 

Pitchline velocity, V, is determined by; 

ir d n p 

V = (18) 

12 

where V is in feet per minute, d is in inches, and n p is in 
revolutions per minute. The tangential component of tooth 
load, W, , is computed from; 
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126000 Pwr 



( 19 ) 



W, = 

n p d 

where W, is in pounds-f orce, Pwr is in horsepower, and d and 
n p are as before. Tooth loading per inch of facewidth is 
computed from; 

Tooth Load per Inch = W t /F (20) 

and the unit lead, a normalized value of the load per inch 
above, is; 

W, Pn* 

Unit Load = (21) 

F 

where the unit load is expressed in pounds-force per square 
inch . 

Mesh frequencies provide information on how often a 
tooth is loaded. Mesh frequencies for parallel axis gears 
are determined by; 

N p n p 

f = (22) 

60 

with f expressed in Hertz. For epicyclic gears, the follow- 
ing are used: 
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NP N R 

(a) f 9 = n s 

N r+ N s 

n r n s 

(b) fp — n s 

Npln N r + N s 

NP Ng 

(c) f n = n s 

N„+ Ng 



(d) f 8 = NP n s 



N s 

(e) f P = 2 n s (23) 

Nrln 

N s 

(f) f„ = NP n 

N„ 



where (a) through (c) are for planetary arrangements and (d) 
through (f) are for star arrangements. 

The K-factor is computed for reference purposes by; 



». (n G + 1) 

K (24) 

F d m G 



The contact stresses are computed according to reference 2 
by; 




(25) 



Bending stresses are computed according to reference 3 by; 



W, P d K s K m 

K v F J 

Individual torques, T, are found by; 



(26) 
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T 



W, d 



(27) 



2000 

while the total output torgue is computed by; 

63 S HP 

T = (28) 

n p 

where T has the units of Thousands of inch-pounds-force in 
both cases. Shaft horsepower, SHP, is the total power 
transferred to the output shaft. 

Weight and size estimates are based on empirical rela- 
tions obtained from a limited number of actual designs. The 
relations used are; 

Weight = Cl[£(d2F) ] c 2 
Length = C3 • £ F 

(29) 

Width = C4 • D 
Height = C5 ■ D 

where the constants used are found in Table 1. All dimen- 
sions are in inches and the weight is in pounds-force 
rounded to three significant figures. 
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Table 1: 



Empirical Constants for Weight and 
Size Formulations 



Constant 


Parallel Axis 


Eoicvclic 


C 1 


1196 .0 


0. 905 


C 2 


0. 34 


0.89 


C 3 


2. 26 


2.85 


C 4 


1. 20 


1.30 




1. 37 


— 


C5 


1. 28 


1 . 20 


D 


Bull Gear 


Ring Gear 




Diameter 


Diameter 



first 


C4 : 


for 


single 


power 


inputs 


second 


C4 : 


for 


double 


power 


inputs 



V. COMPUTATIONAL SUBPROGRAMS LIBRARY FORMULATIONS 
The formulations provided below are for the major compu- 
tational subprograms in Module Four. Those that are self- 
explanatory or are nor computational in nature are only 
described in general. 

t 

A. ARCCOS AND ARCSIN 

These function subprograms find the arc cosine and arc 
sine, respectively, for any two arguments. They were added 
for convenience since not all compilers have them as inrer- 



4 1 



nal functions 



B. AGMAE1 



This function subprogram returns the value of the con- 
stants H, L, and M required for the determination of the 
stress concentration factor, K, , according to reference 1, 
for use in computing the strength geometry factor, J. Taole 
E- 1 in reference 1 provides the tabulated data necessary to 
perform a LaGrangian interpolation for each constant for a 
specified normal pressure angle in degrees. The interpola- 
tion formula used is; 

(On-20) (0 n -14.5) (0„-14.5) (0 n -25) 

F(0J = Fl + p 2 

57.75 -27.50 

(0 n -14.5) (0„-2O) 

+ F 3 (30) 

52.50 

where F represents the appropriate values of H, L, or M. 

C. CKDATA 

FUNCTION CKDATA is called by subroutine AGMA to allow 
the user to selectively change the preprogrammed constants 
by checking if the value entered is zero. If it is zero, 
the current value of the specified constant is not changed. 
This provides for flexibility in changing constants with 
multiple values, and it guards against inadvertantly enter- 
ing a value of zero. 
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D. PO WE RB AND POWERH 

These function subprograms are used to compute the maxi- 
mum allowed service power, in horsepower, that can be trans- 
mitted by a gear according to references 2 and 3. The 
formulation based on the strength rating is; 



n d K v F J s ac K l 

P = 

126000 SF K 0 K m K s P d K„ K T 



(31) 



and the durability raring formulation is; 



n d 



I C v 



126000 SF C, C, C 0 C, 






C, C, 



C R C T 



(32) 



where J and I are the respective geometry factors, F is the 
effective facewidth, n is the speed of d in revolutions per 
minute, and d is the pinion pitch diameter for parallel axis 
or is the sun pitch diameter for epicyclics. All other val- 
ues are the preprogrammed constants. 



E. RTF NDR AND FALFA 

The function subprogram RTFNDR , a slightly modified ver- 
sion of FUNCTION ZEROIN [Ref. 6], is used to find the value 
of the root of the equation programmed in function FALFA. 
This root is required by the subroutine GFJ for the computa- 
tion of the strength geometry factor, J. 
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F. SHRLD 



This function subprogram computes the load sharing ratio 
used in computing the geometry factors. The load sharing 
ratio, m N , is determined by; 

p N tt cos O n 

m N = = (33) 

. 95 Z .95 Z ?nd 

where Z is the length of action defined as; 

z =5 -(V D o 2 -0?" -Vd 2 -^ -Vd 2 -^) (3“) 

The subscripts on the pitch diameters are: 

(1) o : outside diameter; d 0 = d + (2/P d ) 

(2) b : base diameter; d b = d cos 0 t 

For epicyclics, replace the outside diameters is equation 34 
with inside diameters ; d, = d - (2/P d ) . 

G. THICK 

FUNCTION THICK returns the value of the normal arc 
thickness of a tooth at a specified diameter given a thick- 
ness at another diameter. For external gears; 

t 2 = d 2 ( (t, /d, ) + inv 0, - inv 0 2 ) (35) 

and for internal gears; 

h 2 = d 2 ( (t, /d, ) - inv O, + inv 0 2 ) (36) 
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where the subscript 2 represents the desired point and 
subscript 1 represents the known point. The involute func- 
tion is defined as; 



inv x = tan x - x (37) 

The arguments of the involute functions in eguations 35 and 
36 are the transverse pressure angles at the points under 
consideration. The pressure angle at the desired point is 
defined as; 

d, cos 0, 

cos 0 2 = (38) 

The known point is usually taken at the pitch circle where 
d, = d, 0 1 =0 n , and t, is defined as 

P n 7T 

t = = cos Ip (39) 

2 2 P d 



H. GFI 

This subroutine is used to compute the AGMA durability 
geometry factor, I, in accordance with reference 2. The 
geometry factor is defined as; 



I 



cos 0, sin 0 t 



2 



® G 

(n» G ±1) 



( 40 ) 
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where tn N is computed by function SHHLD described above. The 
plus sign applies to external gears and the minus sign 
applies to internal gears. 

I. GFJ 

SUBROUTINE GFJ is used to compute the AG a A strength 
geometry factor, J, in accordance with reference 1 with one 
major difference: the values used are from analytical devel- 
opments and are not scaled to a normal diametral pitch of 
one as are the values used in a graphical layout discussed 
in reference 1. The strength geometry factor is defined as; 



Y c cos 2 ib 

j = (41) 

K, m N 



The load sharing ratio, m N , is computed in FUNCTION 
SHRLD. The stress concentration factor, K, , is determined 
f rom ; 



K, 




(42) 



where H, L, and M are determined in FUNCTION AGHAE1. The 
value of the root fillet radius, r, , is; 



(b - r T ) 2 

(d/2cos 2 ip ) + (b - r T ) 



( 43 ) 



with the dedendum, b = 1.25/P d , and the root tip radius. 
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r T so.28/P nd . The values of t and h are determined from the 
anlytical geometry of the tooth form layout described below. 
The tooth form factor, Y , is defined as; 



Yc 



Pnd 



cos 0 Ln 
cos 0„ 





(44) 



where t and x are also from the tooth form layout mentioned 
previously. The helical factor, C h , is defined as; 



C h 



1 - 




(45) 



where tan v = tan ip sin 0 n for ip < 50°. 
pressure angle at the tip of the tooth, 
figures 6 and 7 and is given by; 



The normal load 
0 Ln , can be seen in 



<Z> 



Ln 




± 



do 



(46) 



where the subscript o pertains to the point on the outside 
diameter and subscript b pertains to the base circle. The 
plus sign applies to internal gears and the minus, to exter- 
nal. The thickness, t 0 , at the outside diameter is deter- 
mined by function THICK. For internal gears, replace the 
outside values with the inside values as before. 

The graphical tooth form layout is a method by which the 
variables h, t, and x can be determined from actual 
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measurements of a tooth form drawn and scaled for a normal 
diametral pitch of one for the case where tooth loading is 
at the tip. Loading at the tip of the tooth is the general 
practice for considering loads on helical gears. Refer to 
Figure 7 for the meanings of h, x, and T where t = 2T. 

Before determining h , x, and T analytically, several refer- 
ence parameters must be determined as suggested by Mclntire 
and Lyon [Ref. 9]. The first is the radius from the center 
of the gear to the tip of the inscribed Lewis stress para- 
bola which is point E in Figure 7. This point is the inter- 
section of the line of action of the tip load, tangent to 
the base circle, with the tooth centerline. The radius to 
this point is; 

dy d b 

r v = = (47) 

2 2cos 0 Ln 

An additional reference point is required to fix the geom- 
etry. The center of of the root fillet is taken as this 
point which can be obtained by a very close approximation. 

To locate this point, the gear center is taken as the origin 
of a cartesian coordinate system with the tooth centerline 
as the vertical axis. Two possible cases exist for the 
location of this point with respect to the base circle. 
Figure 8 shows the case where the point is inside the base 
circle and Figure 9 shows the case where it is outside. The 
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coordinates of this point, (XC,YC), can be found from Fig- 
ures 3 and 9. For both cases it can be seen in Figures 8 
and 9 that; 



HYP = d R + r, (48) 



where d R 


= d 


- 2b = d - 


(2.5 /P d ) . From 


Figure 6, the angle. 


e , is; 




e = inv 0 


t c 

+ sin -1 

d 


(49) 


where t e 


is t 


he chordal 


tooth thickness 


given by; 






T* — ^ « 


t 3 COS 2 i p 


(50) 






U c — W 


6 d 3 


and t is 


the 


normal arc 


tooth thickness 


defined earlier. 


For the 


case 


in Figure 8 


♦ 








XX = 


(HYP) sin e 








xc = 


XX + r, 


(a) 

(51) 

(b) 






YC = ' 


\J HYP 2 - XC 2 



and for the case in Figure 9; 
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(d b /2) 



0, = cos - 1 -- — - 

HIP 

OPP, = (HYP) sin 0, 

OP? 2 s OPP, - r, 

HYP, = V 0PP2 + (d b /2) 2 

(d„/2) 

0 2 = cos ~ 1 

HYP, 

X = 0, ± inv 0 2 - 0 2 

for external gears 

,, + " for internal gears (see Figure 11) 
6 = \ + € 



XC 


= (HYP) 


sin 6 


(a) 








(52) 


YC 


= (HYP) 


cos 6 


(b) 


With the reference 


values 


of r, , XC, 


and YC determined. 



the values of h, t=2T, and x can be analytically determined. 



From Figure 7; 

XT = r, cos a (53) 

YH = r, sin a (54) 

h = r v - YC + YH = r v - YC + r, since (55) 

T = (t/2) = XC - XT = XC - z, cos a (56) 
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T 



( 57 ) 



Y K = 

tan a 

where a must be determined such than; 



YK = 2h or YK - 2h = 0 (58) 

Substituting equations 53 through 57 into 58 yields; 



F(a) = XC 


- r, cos a - 2tan a (r v -YC + r, sin a) - 0 (59) 


Equation 59 


is the function in FALFA called by RTFNDR to 



solve for a. Once a is determined, h can be determined from 
equation 55 and T and t from aquation 56. To obtain x. 



observe the 


following ; 


and 


y = tan - 1 (h/T) 

% = (V2) - y 

i = T tan T, (60) 



While not precise, the identical methodology is used for 
internal gears. Figures 10 and 11 apply. The expressions 
for infernal gears are given without further development; 

5 1 



h = -r v + YC + r, sin a 
T = XC - r, cos a 
t = 2T 

( 61 ) 

y = tan-* (h/T) 

7,= (tt/2) - Y 

x = T tan 7, 

The values for h r t, and x are now used to determine the 
stress concentration factor, equation 42, and the toorh form 
factor, equation 44, required to compute the strength geom- 
etry factor, J, in equation 41. 
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APPENDIX 3 



LIST OF PARAMETERS 



While it is not practical to list all variables used in 
the formulations or the program, it is useful to provide a 
list of the major variables with a cross-reference between 
the analytical names and the FORTRAN names. A detailed 
listing of each common block is also useful when studying 
the program. 



PARAMETER C ROSS- REFER ENC] 



This 


section 


provides a listing of para 


heir analytical 


and FORTRAN names. 


Math 


FORTRAN 


7 ariable 


Symbol 


Name 


Definition 




AKL 


life factor 


K m 


ARM 


load distribution factor 


K 0 


AKO 


overload factor 


Kr 


AKR 


reliability factor 


K s 


AKS 


size factor 


Kt 


ART 


temperature factor 


K v 


AR7 


dynamic factor 


SF 


SF3 


service factor 


C 


CDE 


center distance (theoret 
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CDP 


(E=epicyclic, P=parallel axis) 


C, CF 


surface finish factor 


C H CH 


hardness factor 


C L CL 


life factor 


Cm CM 


load distribution factor 


Co CO 


overload factor 


C p CP 


elastic properties factor 


C R CH 


reliability factor 


C s CS 


size factor 


C T CT 


temperature factor 


C v CV 


dynamic factor 


SF SFH 


service factor 


ip DHELIX 


helix angle (deg) 


HELIX 


helix angle (rad) 


0, DPHI 


transverse pressure angle (deg) 


PHI 


transverse pressure angle (rad) 


0„ DPHIN 


normal pressure angle (deg) 


PHIN 


normal pressure angle (rad) 


D DG 


diameter of gear (in) 


d DP 


diameter of pinion (in) 


d PLN DPLM 


diameter of planet gears (in) 


d„ DR 


diameter of ring gear (in) 



roo-c diameter of a gear (in) 



d s DS 


diameter of sun gear (in) 


F FACES 


facewidth (in) 


FACEP 


(E=epicyclic , P=parallel axis) 
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F/d 


FBYDE 


f/d ratio ( facewidth/diameter) 




FBYDP 


(E=epicyclic, P=parallel 


axis) 


I 


GEOMI 


duribility geometry factor (pinio 




GI 


durability geometry factor (sun) 


J 


GEOMJG 


strength geometry factor 


(gear) 




GEO MJ P 


strength geometry factor 


(pinion) 




GJS 


strength geometry factor 


(sun) 




GJPL 


strength geometry factor 


(planet) 


K 


KFCTRE 


computed !c-factor 






KFCTRP 


(E=epicyclic, P=parallel 


axis) 


f 


MFE 


mesh frequency (Hz) 






MFP 


(E=epicyclic , P=parallel 


axis) 


Mo 


MGOE 


overall reduction ratio 






MGOP 


(E=epicyclic, P=parallel 


axis) 


»G 


MGE 


stage reduction ratio 






MG? 


(E=epicyclic, P=parallel 


axis) 


s G 


NG 


number of teeth, gear 






HP 


number of teeth, pinion 




NP 


NPLNT 


number of planet gears in 


epicycl 


^PLN 


HPLN 


number of teeth, planet 




»R 


HR 


number of teeth, ring 




Us 


HS 


number of teeth, sun 




D 

-d 


PD 


transverse diametral pitch 


Pod 


PHD 


normal diametral pitch 




V 


PLVE 


pitch line velocity (fpm) 






PL VP 


(E=epicyclic , P=parallel 


axis) 



set 
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Pwr 


PWRE 


power split per gear pair (hp) 




PWRP 


(E=epicyclic, P=parallel 


axis) 


n in 


RPMIN 


source speed input (rpm) 




tout 


RPMOUT 


output shaft/propeller speed (rpm) 


n* 


RPMI 


stage input speed, epicyclic (rpm) 




RPMO 


• stage output speed, epicyclic (rpm) 


n PLN 


RPMPL 


planet speed, epicyclic 


(rpm) 


n P > n G 


RPMP 


stage pinion and gear sp 


eed, parallel 






axis (rpm) 




Sac 


SAC 


allowable contact stress 


number 


Sat 


SAT 


allowable bending stress 


number 


SHP 


SHP 


shaft horsepower, output 


(hp) 


s, 


SIGBE 


bending stress (psi) 






SIG3P 


(E=epicyclic, P=parallel 


axis) 


s c 


SIGHE 


contact stress (psi) 






SIGHP 


(E=epicyclic, ?=paraliel 


axis) 


T 


TORQE 


torgue (k in-lb) 






TORQP 


(E=epicyclic , P=parallel 


axis) 


w« 


WTE 


tangential tooth load (lb) 



WTP (E=epicyclic, P=parallel axis) 



II. COMMON BLOCK DETAILS 

The following provides information concerning the vari- 
ables in each common block. The numbers in parentheses are 
the size of the array where applicable. 
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COMMON 


3L0CK AGMA3 


(FOR STRENGTH RATING) 


SFB 


• 

• 


R*4 (2,2) ; 


service factor 


AKV 


• 

« 


R* 4 ; 


dynamic factor 


AKS 


« 


R*4 ; 


size factor 


AKM 


« 

• 


R*4 ; 


load distribution factor 


AKO 


• 

• 


R*4 (2) ; 


overload factor 


SAT 


« 

# 


H*4 (6) ; 


allowable bending stress number 


AKL 


• 

• 


R*4 (2) ; 


life factor 


AKR 


• 

# 


H^4 (6) ; 


reliability factor 


AKT 


• 

• 


R a 4 ; 


temperature factor 



COMMON 


BLOCK 


AGMAH 


(FOR DURABILITY RATING) 


SFH 


: R*4 


(2,2) ; 


service factor 


cv 


: R*4 


(3) ; 


dynamic factor 


cs 


: R*4; 




size factor 


CM 


: R*4 


(2) ; 


load distribution factor 


CF 


: R*4; 




surface finish factor 


CO 


: R*4 


(2) ; 


overload factor 


SAC 


: 3* 4 


(6) ; 


allowable contact stress number 


C? 


: R*4 ; 




elastic properties factor 


CL 


: R* 4 


(2) ; 


life factor 


CH 


: R ;i 4; 




hardness factor 


CT 


: R*4; 




temperature factor 


CR 


: R*4 


(6) ; 


reliability factor 
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COMMON BLOCK DESDAT (DESIGN P ARAM ETERS , INPUT) 



PWRIN s 


R*4 


(2) ; 


source power input (hp) 


RPMIN : 


R*4 


(2) ; 


source speed input (rpra) 


RPMOUT: 


R*4 ; 




output shaft/propeller speed (rpm) 


DHELIX : 


R*4 


( 3 ) ; 


helix angle (deg) 


HELIX : 


R* 4 


( 3 ) ; 


helix angle (rad) 


PD : 


R* 4 


( 3 ) ; 


transverse diametral pitch 


PND 


R* 4 


( 3 ) ; 


normal diametral pitch 


DPHI : 


R*4 


( 3 ) ; 


transverse pressure angle (deg) 


PHI : 


R*4 


( 3 ) ; 


transverse pressure angle (rad) 


DPHI N : 


3*4 


( 3 ) ; 


normal pressure angle (deg) 


PHIN : 


R* 4 


( 3 ) ; 


normal pressure angle (rad) 


NDI?P : 


1*4 ; 




number of different power sources 


IARR 


1*4; 




arrangement code (1=parallel axis, 
2=epicyclic) 


I EPIC : 


1*4 


( 3 ) ; 


epicyclic code ( 1 =planetary , 2=star) 


I HARD : 


1*4 


( 3 , 2 ) ; 


hardness range code (1-6, see SU3R. 
AGMA) 


iopro : 


1*4 ; 




operational profile code (1=naval pro 
file full power 5^ max; 2=other, max 
power continuous) 


NPWRIN: 


1*4 ; 




number of power sources (inputs) 


IPWRSR: 


1*4 


(2) ; 


power source code (1=turbine or motor 
2=multicylinder internal combustion 
engine) 


NRED : 


1*4; 




number of reduction stages 
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NPATH 

NPLNT 

NHELX 

COMMON 

HGOE 

MGE 

RPMI 

RPMPL 

RPMO 

PWRS 

DS 

DPLN 

DR 

FACEE 

GI 

GJS 

GJPL 

NS 

NPLN 

NR 

COMMON 

PWRFAC 

MGOP 

MGP 



: 1*4; number of power paths { 1=single , 2=dual) 

: 1*4 (3) ; number of planet gears in epicyclic set 

: 1*4; number of helicies (1 = single, 2 = double) 



BLOCK DSSEPC (EPICYCLIC DESIGN PARAMETERS) 



R*4 ; 




overall reduction ratio 


R*4 


(3) ; 


stage reduction ratio 


R* 4 


(3) ; 


stage input speed (rpm) 


3*4 


(3) ; 


planet speed (rpm) 


R* 4 


(3) ; 


stage output speed (rpm) 


R* 4 


(3) ; 


stage power split per planet (hp) 


R* 4 


(3) ; 


diameter of sun gear (in) 


R* 4 


(3) ; 


diameter of planet gears (in) 


R* 4 


(3) ; 


diameter of ring gear (in) 


R* 4 


(3) ; 


facewidth (in) 


R*4 


(3) ; 


durability geometry factor (sun/planet) 


R*4 


(3) ; 


strength geometry factor (sun) 


R* 4 


(3) ; 


strength geometry factor (planet) 


1*4 


(3) ; 


number of teeth, sun 


1*4 


(3) ; 


number of teeth, planet 


1*4 


(3) ; 


number of teeth, ring 



BLOCK DSSPP.L ( PARALLEL AXIS DESIGN PARAMETERS ) 
: R*4 (2,3) ; stage power split factor 
: R*4 (2) ; overall reduction ratio 
: R*4 (3,2); stage reduction ratio 
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RPilP 



PWRP 

D? 

DG 

FACS? 
GEO MI 
GEOMJG 
GEO M J P 
NP 
NG 

COMMON 
PL VE 
FBYDE 
CDS 
tfTE 
TLPIE 
UNTLDE 
MFE 

KFCTRE 

SIGHE 

SIGBE 

TORQE 

RPME 

PDIAME 

WGHTE 



R'*4 


(6,2) ; 


5*4 


(6 # 2) ; 


R*4 


(3,2) ; 


5*4 


(3,2) ; 


R*4 


(3,2) ; 


R*4 


(3,2) ; 


R*4 


(3,2) ; 


R*4 


(3,2) ; 


1*4 


(3,2) ; 


1*4 


(3,2) ; 



stage pinion and gear speed (rpm) 
srage power split per gear (hp) 
diameter of pinion (in) 
diameter of gear (in) 
facewidth (in) 
duribility geometry factor 
strength geometry factor (gear) 
strength geometry factor (pinion) 
number of teeth# pinion 
number of teeth# gear 



BLOCK RESEPC 



5*4 


(3) ; 


R*4 


(3) ; 


5*4 


(3) ; 


3*4 


(3) ; 


R*4 


(3) ; 


R^4 


(3) ; 


3*4 


(3,3) 


5*4 


(3) ; 


R*4 


(3) ; 


5*4 


(3) ; 


R*4 


(3,3) 


R* 4 


(3,3) 


3* 4 


(3,3) 


R* 4 ; 





f EP ICY CL IC PARAMETERS . RESULTS ) 
pitch line velocity (fpm) 
f/d ratio (facewidth/sun diameter) 
center distance (theoretical) (in) 
tangential tooth load (lb) 
tooth load per in (lb/in) 
unit load (psi) 
mesh frequency (Hz) 
computed k-factor 
contact stress (psi) 
bending stress (psi) 
torgue (k in-lb) 
gear speeds (rpm) 
pitch diameters (in) 
gear set weight estimate (lb) 
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SPCWTS: 


R*4 ; 




specific weight (lb/hp) 


MTHE : 


1*4 


( 3 , 3 ) ; 


tooth numbers 


ISIZES: 


1*4 


( 3 ) ; 


length, width, height estimates (in) 


COMMON BLOCK 


RESPRL 


(PARALLEL AXIS PARAMETERS. RESULTS) 


PLVP : 


3* 4 


( 3 , 2 ) ; 


pitch line velocity (fpm) 


FBYD? : 


R* 4 


( 3 , 2 ) ; 


f/d ratio (facewiath/pinion diameter 


CDP : 


R*4 


( 3 , 2 ) ; 


center distance (theoretical) (in) 


WTP : 


3* 4 


( 6 , 2 ) ; 


tangential tooth load (lb) 


TLPIP : 


R*4 


( 6 , 2 ) ; 


tooth load per inch (lb/in) 


UNTLDP: 


R* 4 


( 6 , 2 ) ; 


unit load (psi) 


MFP : 


R* 4 


( 3 , 2 ) ; 


mesh frequency (Hz) 


KFCT3P : 


R* 4 


( 6 , 2 ) ; 


computed k-factor 


SIGHP : 


R* 4 


( 3 , 2 ) ; 


contact stress (psi) 


SIG3P : 


R* 4 


( 6 , 2 ) ; 


bending stress (psi) 


TORQP : 


3*4 


( 6 , 2 ) ; 


torque (k in-lb) 


PDIAMP : 


3*4 


( 6 , 2 ) ; 


pitch diameters (in) 


SCDMIN; 


R* 4 ; 




minimum source center distance (in) 


SCDMAX: 


R* 4 ; 




maximum source center distance (in) 


SHP : 


3*4 ; 




shaft horsepower, output (hp) 


WGHTP : 


R*4 ; 




gear set weight estimate (lb) 


SPCWTP: 


R*4 ; 




specific weight (lb/hp) 


TRQOUT: 


R*4 ; 




torque, output (k in-lb) 


MTHP : 


1*4 


( 6 , 2 ) ; 


tooth numbers 


ISIZEP: 


1*4 


( 3 ) ; 


length, width, height estimates (in) 
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APPENDIX C 



5 EG AD SAMPLE RUNS 



This appendix contains samples of actual terminal ses- 
sions using HEGAD. For the sake of brevity, only two com- 
plete sessions are included. However, a number of analysis 
and design runs were made using a full range of options and 
configurations, and they compared favorably to actual 
designs. The comparisons are not shown here due to the pro- 
prietary nature of the designs used for verification. The 
first example is an analysis run for a locked train, double 
reduction gear set with two different inputs. Following it, 
is the results section from a design run using the identical 
parameters as the analysis run. The second example is a 
double reduction epicyciic gear set with the complete analy- 
sis session followed by the results section of a design run 
as before. The analysis and design sessions are identical 
with one exception. A seed for a random number generator is 
requested in the design option instead of diameters and 
facewidths as in the analysis option. For those cases where 
an infeasible design is generated, a message will alert the 
user and the program will continue. To obtain a feasible 
design, or just a different one, rerun the program and 
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provide a different seed for the random number generator 
This method was used on several occasions to obtain the 



desired results, 
user can then use 
that more closely 



Once a feasible design is obtained, the 
the analysis option to obtain a design 
suits his needs. 
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CHOOSE OPERATIONAL PROFILE CODE BELOW: 
OPERATIONAL MODE SERVICE PROFILE CODE 
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POWER SPLIT HP | 21250.1 10625 . | 10625. | 21250 
SPEED RPM | 5980. J 1851. | 1851. | 300 
NUMBER OF TEETH | 76 | 247 | 73 | 452 
NORMAL DIAMETRAL PITCH I 5.493 | 4.273 
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CHOOSE OPERATIONAL PROFILE CODE BELOW: 
OPERATIONAL MODE SERVICE PROFILE CODE 
FULL POWER 5 PERCENT MAX 1 

MAXIMUM LOAD CONTINUOUS 2 



•'«* ENTER OPERATIONAL FROFILE CODE: 
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WHICH DIAMETRAL PITCH WILL YOU SPECIFY? 
(1=TRANSVERSE, 2=N0RMAL) : 
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ENTER HARDNESS CODES FOR SUN/PLANETS AND RING (HCSUN,HCRING) : 
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ENTER HARDNESS CODES FOR SUN/PLANETS AND RING (HCS UN , HCRI NG ) : 
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FUNCTION AGMAE1 (PHIN,II) 4MOD0010 

4MOD0020 

CODED BY: LT J.L. PAQUETTE, USN JAN 1982 4MOD0030 

NAVAL POSTGRADUATE SCHOOL MONTEREY, CA 93940 4MOD0040 

4MOD0050 
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APPENDIX E 



FIGURES 








Figure 1: Single Reduction Parallel Axis Arrangements 
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Double Power Path 

Figure 2: Double Reduction Parallel Axis Arrangements 
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Figure 3: Triple Reduction Parallel Axis Arrangements 
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Figure 4: Single Reduction Epicvclic Arrangements 

(from Ref. 7) 



176 




Start 




Figure 5: Flow Chart of the 3EGAD Package 
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Figure 6: External Tooth Dimensions 
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Figure 9: Fillet Center Location - Outside Base Circle 
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Figure 10: Internal Tooth Dimensions 
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